Dualities in population genetics: a fresh look with 

new dualities. 



Gioia Carinci'"', Cristian Giardina'''' , 
Claudio Giberti*''^ Frank Redig<">. 

^'^^ Department of Mathematics, University of Modcna and Reggie Emilia 
via G. Campi 213/b, 41125 Modena, Italy 
Department of Engineering, University of Modena and Reggio Emilia 
via Giovanni Amendola 2, 42122 Reggio Emilia, Italy 
('^^ Delft Institute of Applied Mathematics, Technische Universiteit Delft 
Mekelweg 4, 2628 CD Delft, The Netherlands 

February 14, 2013 



Abstract 

We apply our general method of duality, introduced in [6] , to models 
of population dynamics. The classical dualities between forward and 
ancestral processes can be viewed as a change of representation in the 
classical creation and annihilation operators, both for diffusions dual 
to coalescents of Kingman's type, as well as for models with finite 
population size. 

Next, using SU{1, 1) raising and lowering operators, we find new 
dualities between the Wright-Fisher diffusion with d types and the 
Moran model, both in presence and absence of mutations. These new 
dualities relates two forward evolutions. From our general scheme we 
also identify self-duality of the Moran model. 
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1 Introduction 



In the works [6l [7] duality between two stochastic processes, in the context 
of interacting particle systems and non-equilibrium statistical mechanics, 
has been related to a change of representation of an underlying algebra. 
More precisely, if the generator of a Markov process is built from lowering 
and raising operators (in physics language creation and annihilation oper- 
ators) associated to a Lie algebra, then different representations of these 
operators give rise to processes related to each other by duality. The inter- 
twiner between the different representations is exactly the duality function. 
Furthermore, self-dualities can be found using symmetries related to the 
underlying Lie algebra (see also [15]). 

The fact that generators can be built from raising and lowering oper- 
ators is a quite natural assumption. E.g., in interacting particle systems 
dynamics consists of removing particles at certain places and putting them 
at other places. If the rates of these transitions are appropriately chosen, 
then the operators of which the effect is to remove or to add a particle (with 
appropriate coefficients) form a Lie algebra. For diffusion processes the gen- 
erator is built from a combination of multiplication operators and (partial) 
derivatives. For specific choices, these correspond to differential operator 
representations of a Lie algebra. If this Lie algebra then also possesses a 
discrete representation, then this can lead to a duality between a diffusion 
process and a process of jump type, such as the well-known duality between 
the Wright-Fisher diffusion and the Kingman's coalescent. 

It is the aim of this paper to show that this scheme of finding dual 
processes via a change of representation can be applied in the context of 
mathematical population genetics. First, we give a fresh look at the clas- 
sical dualities between processes of Wright-Fisher type and their dual co- 
alescents. These dualities correspond to a change of representation in the 
creation and annihilation operators generating the Heisenberg algebra. For 
population models in the diffusion limit (infinite population size limit) the 
duality comes from the standard Heisenberg generators x, d/dx and another 
discrete representation, known as the Doi-Peliti representation. The inter- 
twiner is in this case simply the function D{x,n) = x". In the case of finite 
population size, dualities arise from going from a finite-dimensional repre- 
sentation (finite dimensional creation and annihilation operators satisfying 
the canonical commutation relations) towards the Doi-Peliti representation. 
The intertwiner is exactly the hypergeometric polynomial found e.g. in [lOj . 
[9], and gives duality between the Moran model with finite population size 
and the Kingman's coalescent. 

Next we use the SU{1, 1) algebra to find previously unrevealed duali- 
ties between the discrete Moran model with (or without) mutation and the 
Wright Fisher diffusion with (or without) mutation, as well as self-duality 
of the discrete Moran model with mutation. 
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These are in fact applications of the previously found dualities between 
the Brownian energy process and the symmetric inclusion process, as well as 
the self-duality of the symmetric inclusion process, which we have studied in 
another context in [6], [7], [8]. Put into the context of population dynamics, 
these dualities give new results for the multi-type Moran model with or 
without mutation, as well as the multi-type Wright Fisher model with or 
without mutations. 

The rest of our paper is organized as follows. In section [2] we give the 
general setting and view on duality. Though many elements of this formal- 
ism are already present in previous papers, we find it useful to put these 
together here in a unifying, more transparent and widely applicable frame- 
work. In section [3] we discuss dualities in the context of the Heisenberg 
algebra. This leads to dualities between diffusions and discrete processes, 
as well as between different diffusion processes, and finally between different 
discrete processes. In section [J] we show how this gives the dualities between 
forward (in time) population processes and their ancestral dual coalescents. 
In section[5]we apply the SU{1, 1) algebra techniques, developed previously 
in [7] in the context of population dynamics, finding new dualities, this time 
between two forward (in time) processes: duality between the Wright-Fisher 
diffusion and the Moran model, and duality of the Moran model with itself. 
We give three concrete computations using these new dualities as an illus- 
tration. 



2 Abstract Setting 

2.1 Functions and operators. 

Let O, O be metric spaces. ^(O), '^(0), "^($7 x Q) denote the sets of bounded 
continuous real- valued functions on fi, fi, x Jl. In what follows, this choice 
of function space is not stringent, and can be replaced if necessary by other 
function spaces such as spaces or Sobolev spaces. 

For a function x/j : Q x Q M. and linear operators K : T){K) C ^($7) — )• 
"^{n), k : Ti{k\ C ^(O) ^ M we define the left action of K on V and the 
right action of X on -0 via 

(K,V)(x, y) = {Kij{;y)) (x) , {Kri^){x, y) = {k^{x, ■)){y) , (1) 

where we assume that -0 is such that these expressions are well-defined, i.e. 
V'(-,y) eD(K),V'(x,-) e^(^)- 

An important special case to keep in mind is O = {l,...,ri} and $7 = 
{1, . . . , m} finite sets, in which case we identify functions with column vec- 
tors. In that case, an operator K (resp. K) on functions on (resp. Jl) 
coincides with a n x n (resp. m x m) matrix. Denoting by K (resp. K) 
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such matrices, one has 



Kiilj{x,y) = ^^{x,z)'4){z,y) = (IK^/;)(x,y) , 



(2) 



(3) 



Namely, left action of K corresponds to left matrix multiplication and right 
action of K corresponds to right multiplication with the transposed matrix. 
The same picture arises when 0, are countable sets. 

For two operators Ki,K2 working on the same domain we denote, as usual, 
K1K2 their product or composition, i.e. 



We further abbreviate =^(0) for algebras of linear operators working on a 
common domain of functions in ^(Q). These algebras are often representa- 
tion of an abstract algebra (such as e.g. the Heisenberg algebra). For a 
general algebra M', we define the dual algebra M'* as the algebra with the 
same elements as in M' but with product "*" defined by a * 6 = 6 • a, where 
• is the product in M". A typical example in the finite dimensional setting 
is the algebra of n x n matrices, where the map maps the algebra 

into the dual algebra {[ABY = B'^A^). 

Definition 2.1. For a function il) : n x ^ "M. we say that ip is left 
exhaustive if the relation Kitp = implies -fC = 0, and correspondingly we 
call ip right exhaustive if the relation K^ip = implies K = 0. Notice that 
in the context of finite sets 0,, Q this just means that the matrix associated 
to ip is invertible. 

2.2 Duality 

Definition 2.2. Let K G .^{9), K e .^{fl) and D : n x n ^ R. Then 
we say that K and K are dual to each other with duality function D if 



We denote this property by K — )■ K. 

In the following we collect elementary but important properties of the rela- 
tion — 7>^. The proof of these properties is elementary and left to the reader. 



Theorem 2.1. Let Ki,K2 G =^(0), Ki,K2 e ^{h). Suppose that Ki 
Ki , K2 — K2 , and further ci , C2 G M then we have 

L ki -^^ Ki, with D{x,y) = D{y,x). 



{KiK2f){x) = {K^{K2f)){x) . 



KiD = KrD . 



(4) 



4 



2. ciKi + C2K2 ciKi + C2K2. 

3. K1K2 -^^ K2K1, in particular i^f -^^ , n e N. 

4- [Ki,K2] —^^ [K2,Ki] = —[Ki,K2], i.e., commutators of dual opera- 
tors have opposite sign. 

5. If S e ^{Q,) commutes with Ki, then Ki -^^'^ Ki, and if S e =^(0) 
commutes with Ki, then Ki -^'^^D 

6. If for a collection {Ki,i G /} C ^{Q.), and {Ki : i e 1} C .^{Q) 
we have Ki — 7>^ Ki then every element of the algebra generated by 
{Ki,i € /} is dual to an element of the algebra generated by {Ki : i G 
/}. More precisely: 

ii • • • ifc ik ' ' ' h 

for all ni, . . . ,ni^ G N, and for constants {q : i G /} 

Y^cK^^^'Y^aK, . 

7. Suppose ^only in this item) that Ki G ^(f^i), Ki G K2 G 
^{^2), K2 e^^{n2). If Ki ^-^1 Ki and K2 K2, then Ki O 
K2 ^■^i®^^ K10K2, where Di0D2{xi,X2;yi,y2) = Di{xi,yi)D2{x2,y2). 

8. If K and K generate semigroups St = e*^ and St = e*^ then K — K 
implies St — St. If moreover these semigroups correspond to Markov 
processes {Xt, t > 0}, {Xt, t > 0} on Q, Q, the relation St — St reads 
in terms of these processes: 

E,D{Xt,x)=EsD{x,Xt) , (6) 

for all X ^ Vt,x ^ n and t > 0. On the left-hand side of the above 
equation we have expectation with respect to the the law of the {Xt}t>o 
process started at x, while on the right-hand side we have expectation 
with respect to the law of the {Xt}t>o process initialized at y. 

Remark 2.1. Item 6 of theorem \2.1\ is useful in particular if the collection 
{Ki,i G /} C c'^{il.) is a generating set for the algebra. Indeed, then every 
element of the algebra has a dual by ([5]), and it suffices to know dual opera- 
tors for the generating set to infer dual operator for a general element of the 
algebra. In practice, one starts from such a generating set and the commu- 
tation relations between its elements (defining the algebra) and associates to 
it by a single duality function a set of dual operators with the same commu- 
tation relations up to a change of sign. In other words, one moves via the 
duality function from a representation of the algebra to a representation of 
the dual algebra. 
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(5) 



Remark 2.2. The relation ([6]) is the form in which one usually formulates 
duality between two Markov processes. Remark however that the relation 
St St between the semigroups is more general. It may happen that St 
is a Markov semigroup, whereas St is not. E.g., mass can get lost in the 
evolution according to the dual semigroup St, which means Stl ^ 1, or it 
can happen that St is not a positive operator. 

In [13], the author studies the so-called duality space associated to two 
operators. In our notation, this is the set 

D{K, K) = {D -.nxn^R: K K} . (7) 

Our point of view here is to consider for a fixed duality function the set of 
pairs of operators (K, K) such that K — K. These operators then usually 
form a representation and a dual representation of a given algebra, with D 
as inter- twiner. 

In theorem 12.11 item 5, we see that we can produce new duality functions 
via "symmetries", i.e., operators commuting with K or K. In the context 
of finite sets $7 = $7 = {1, . . . ,n}, we have more: if there exists an invert- 
ible duality function D, then all other duality functions are obtained via 
symmetries acting on D. We formulate this more precisely in the following 
proposition. We will give examples in the subsequent sections. 

Proposition 2.1. 1. Let VL = Q. = {1, . . . ,n}, and let K -^^ K. Sup- 
pose furthermore that the associated nxn matrix D is invertible. Then, 
if K — K , we have that there exists S commuting with K such that 
D' = SD. 

2. For general $7, $7' we have the following. Suppose D and D' are duality 
functions for the duality between K and K . Suppose furthermore that 
D' = SiD, for some operator S. Then we have 

{KS - SK)iD = . (8) 

In particular if D is left exhaustive, then we conclude [S.,K] = 0. 
Similarly, if D' = SpD then 

{KS - Sk)rD = . (9) 

and if D is right exhaustive, then we conclude [S, K] = Q . 

Proof. In the proof of the first item, with slight abuse of notation, we 
use the notation K, K, D, S both for the operators and for their associated 
matrices. We have, by assumption 

KD = DK'^. 
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By invertibility of D, S = D'D ^ is well-defined and we have 
SK = D'D-^K = D'K'^D-^ , 

and 

KS = KD'D-^ = D'K^D-^ , 

hence [K, S] = 0. 

For the second item, use K K,K K to conclude 
{KS)iD = Ki{SiD) = KiD' = K^D' , 

as well as 

{SK)iD = SiKriD) = KriSiD) = Kr{D') , 

Hence 

{[S,K])iD = 0. 

The remaining part of the proof (for the right action case) is similar. □ 
2.3 Self duality 

Self-duality is duality between an operator and itself, i.e., referring to def- 
inition 12.21 Q = Q and K = K. The corresponding duality function such 
that K —7-^ K is then a function D : x $7 — )• M and we call it a self-duality 
function. For self-duality, of course, all the properties listed in theorem 12.11 
hold. 

In the finite case O = {1, . . . ,n}, a self-duality function is a n x n matrix 
and self-duality reads, in matrix form, 

KL> = D^F . 

Therefore, in this setting such a matrix D can always be found because every 
matrix is similar to its transposed [16], i.e., self-duality always holds with 
an invertible D. 

Other self-duality functions can then be found by acting on a given self- 
duality function with symmetries of K (i.e. operators S commuting with K), 
as we derived in item 5 of theorem 12.11 and in proposition 12.11 In particular 
we have in this finite context, in the notation ([7]): 

D(i^, K) = {SD : [S, = 0} , (10) 

with D an arbitrary invertible self-duality function. So this means that the 
correspondence between self-duality functions and symmetries of K is one- 
to-one. This characterization of the duality space has the advantage that 
the set of operators commuting with a given operator is easier to identify. 
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In the finite setting, if K is tlie generator (resp. transition operator) of 
a continuous-time (resp. discrete-time) Markov chain, then if this Markov 
chain has a reversible probabihty measure : — )• [0, 1], a duahty function 
for self-duahty is given by the diagonal matrix 

This is easily verified from the detailed balance relation //(x)IK.(3;, y) = 
K{y,x)^{y). This "cheap" duality function is usually not very useful since 
it is diagonal, but it can be turned in a more "useful" one by acting with 
symmetries. All the known self-duality functions in discrete interacting par- 
ticle systems such as the exclusion process, independent random walkers, 
the inclusion process, etc. can be obtained by this procedure [6]. 

3 Dualities in the context of the Heisenberg alge- 
bra. 

In this section we focus on representations of Heisenberg algebra and its 
dual algebra, and the corresponding duality functions that connect these 
different representations. 

3.1 Standard creation and annihilation operators 

As a first example, let us start with the operators yl^, ^4 working on smooth 
functions / : R — )■ M with compact support, defined as 

Af{x) = f'ix), A^fix) = xf{x) . (11) 

These operators satisfy the canonical commutation relation [A, A^ = I, 
where I is the identity operator. 

The same commutation relation, up to a negative sign, can be achieved using 
operators working on discrete functions. Considering 

af{n) = nf{n-l), f{n) = f{n + I) , (12) 

acting on functions / : N M, we have [a, = —I. Therefore, in view 
of the item 4 of theorem 12.11 the operators a, are natural candidate for 
duality with ^, 

To find D such that A — >^ a, we use the definition 12.21 

AiD{x, n) = D'{x, n) = arD{x, n) = nD{x, n — 1) , 

which yields 

D{x,n) = ^ rijcn-k 
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with {ci : i € N} a sequence of constants. In the same way, the duahty 
condition A^^ -^^ , produces 

AjD{x, n) = xD{x, n) = a].D[x, n) = D{x, n + 1) , 

which gives 

D{x,n) = x'^D{x,0) , 

with D{x,0) an arbitrary function. Therefore, if we want both duahties to 
hold with the same duahty function, then we are restricted to the choice 

D{x, n) = cox" . 

Without loss of generality we can choose cq = 1. 

As a consequence, by using item 6 of theorem l2.1l we obtain the following 
result. 

Theorem 3.1. For < n < m, /ei a„ : M — )• M 6e a finite sequence of 
polynomials. A differential operator of the form 

m m 
n=0 n=0 

is dual with duality function D{x,n) = to 

m 
n=0 

where the operators a, are defined in H^). 
3.2 Generalization 

In the following proposition we show how to generate the duality functions 
for more general generators A, A^ of a representation of the Heisenberg al- 
gebra, namely by repetitive action of the creation operator on the "vacuum" 
which is annihilated by the operator A. 

Proposition 3.1. Suppose [A, A^] = /, and let D{x,n) he functions such 
that 

{AlrD{x,0)=D{x,n), 

AiD{x,0) = 0, (13) 

then A a and A^ -^^ a^, where a,a^ are the discrete representation 
defined in (fT^ . As a consequence, for a finite sequence of polynomials an, 
with < n < m, we have the analogue of theorem \3.1]: 

m m 

^a„(At)A"^^ J] a"a„(at). 

n=0 n=0 
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Proof. We have — by the assumption on in (|13p and the 
definition of in (jl2p . We therefore have to prove A — t-^ a. Start from the 
commutation relation [A, A^] = / to write 

A(At)" = (At)'^A + [A, (At)"] = (At)"A + n(At)"-i. 

Then use the assumptions (fT3|) to deduce 

A;i^(x,n) = Ai{alrD{x,0) = Ai{AlrD{x,0) 

= (A|)"AiL'(x, 0) + n(A|')"-ii?(x, 0) = n(A|)"-ii?(x, 0) 

= n(aJ)"-^L'(x,0) 

= nD{x, n — 1) = arD{x, n) . 

□ 



As an apphcation, we can choose hnear combinations of multiphcation 
and derivative 

A = cix + C2-^, A''' = C3X + C4-^, (14) 
ax ax 

with the real constants satisfying C2C3 — C1C4 = 1, then we satisfy the com- 
mutation relation [A, A^] = /. To find the corresponding duality function 
that "switches" from A, A^ to a,a\ we start with 

AiD{x, 0) = cixD{x, 0) + C2D'{x, 0) = a^L>(x, 0) = 

which gives as a choice 

Z)(x,0) = expf-^^ 



C2 2 



and next, 



(J \ n 
C32; + C4— j D(x,0). 



An important particular case is when ci = C2 = 1/2 and C3 = — C4 = 1. With 
this choice one finds that the duality function is D(x,n) = ^'^Hn{x), 
where Hn is the Her mite polynomial of order n. 



3.3 Dualities with two continuous variables 

Within the scheme described in section [2] we can also find dualities between 
two operators both working on continuous variables, as the following exam- 
ple shows. 
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Consider again the operators A, in (jlip . A "dual" commutation relation 
(in the sense of item 4 of theorem 12. ip can be obtained by considering a copy 
of those operators and exchanging their role. Namely, we look for dualities 
d/dx — >^ y,x -^^ d/dy. Imposing that the left action of d/dx (resp. x) 
does coincide with the right action of y (resp. d/dy) one finds the duality 
function D{x, y) = e^^ . As a consequence one immediately has the following 

Theorem 3.2. For < n < m, /et a„ : M — ;> M 6e a finite sequence of 
polynomials. A differential operator working on smooth functions of the real 
variable x and with the generic form 

n=0 

is dual, with duality function D{x,y) = e^^ , to the operator working on 
smooth functions of the real variable y given by 

m , 
n=0 ^ 

If we specify that the operators work on functions / : [0, oo) — )• M and 
when the operators can be interpreted as generators of diffusions one has 
the following 

Corollary 3.1. The diffusion generator 

^ = (CIX^ + C2X)^ + {C3X) — 

with ci > 0, C2 > is dual to 

^ = c,y—^^{c,y +c,y)- 

with duality function D{x,y) = e^^. For the corresponding diffusion pro- 
cesses {Xt : t > 0}, {Yt : t > 0} we thus have 

The particular case C2 = gives is that the diffusion generator c\x^^ jdx^ + 
c-ixdjdx is self-dual. 

3.4 Discrete creation and annihilation operators 

The following example starts from a finite dimensional representation of 
the Heisenberg algebra (in the spirit of O [9]). We consider = Ojy = 
{0, . . . , A^} and $7 = N. For functions / : Oat ^ M we define the operators 

a;v/(fe) = (A - fe) /(fc + 1) + (2fe - AT) /(fc) - fc/(fc - 1) , 

r=0 \k) 
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with the convention /(—I) = f{N + 1) = 0. Consider 



n (kn)- - fe(fc-l)...(fc-(n-l)) 



with the convention D^ik, 0) = 1, D]si{k, + 1) = 0. Let us denote by Wn 
the vector space generated by the functions k i— )■ Di\i(k,n),0 <n<N. 

Proposition 3.2. 

{aN)iDj\[{k,n) = nDN{k,n — 1), V 1 < n,V k > n — 1 , 
MiDNik, 0) = OVO<A;<iV, 

(ajv)iL>Ar(A;,n) = DN{k,n + 1) y < n < N,k > n . (17) 
As a consequence, as operators on we have 

[aAr,a]v]=/, (18) 

i.e., aTVjajy form a finite dimensional representation of the canonical com- 
mutation relations. 

Proof. Straightforward computation. □ 

Remark 3.1. Notice that in the limit N ^ oo, putting k/N = x, and 
f{k) = (j){x) = (j){k/N), aNf{k) converges to d(j)/dx. Next, notice that 
Dj\l{k,n) = (/>^^(x), where 



4")(.) = x(^)...f^^ 




converges to x". The effect of the operator ajy on D]\f{k,n) is to raise 
the index n by one. Since Diq{k,n) = cl)^^\x) — )• for N ^ oo and 
a^j^D]^{k,n) = j4^(/)^^(x) — > x'^'^'^ for N — > oo, we conclude that in the 
limit iV — )• oo, the operator ajy coincides with the multiplication operator A'^ 
defined in (|lip . 



The discrete finite dimensional representation of the Heisenberg algebra 
given in proposition 13.21 will be used at the end of section U] to fit within 
the scheme of a change of representation the classical duality between the 
Moran model and the block-counting process of the Kingman coalescence. 
We end this section with a comment on relation between the discrete repre- 
sentation and Binomial distribution. This also offers an alternative simple 
way to see the commutation relation (jlSp . 
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In many models where there is duahty or self-duality, there exists a 
one-parameter family of invariant measures fp, and integrating the dual- 
ity function w.r.t. these measures usually gives a simple expression of the 
parameter p. In the context of diffusion processes with discrete dual, this 
relation is usually that the duality function with n dual particles integrated 
over the distribution Vp equals p^. A similar relation connects the polyno- 
mials Di\f(k,n) in (|16|) to the binomial distribution. This general relation 
between a natural one-parameter family of measures and the duality func- 
tions cannot be a coincidence and requires further investigation. 

The polynomials D^[k,n) are (as a function of k) indeed naturally as- 
sociated to the binomial distribution. Denoting by 

the binomial distribution with success probability p G [0, 1], we have 

AT 

Y,DN{k,n)yM,p{k)=p'" ■ (19) 
For a function / : f^^r ^ M we define its binomial transform 3/^ f : [0, 1] — ?> M 

by 

N 

{^f){p) = Y.f{k)uMAk). (20) 
If for such /, we write its expansion 

N 

f{k) = Y,CrDN{k,r) 

r=0 

we say that / is of degree / if 7^ and all higher coefficients Ck,k > I are 
zero. We then have, using (fT9]l and (pO|) . 

N 
r=0 

The function / and have therefore the same components with respect 
to two different bases: one given by {Di\f{k, r), r = 0, . . . , A'^} which is a base 
of M^"*"^ and the other given by {p^, r = 0, . . . , N} which is a base of the 
space of polynomials on [0, 1] of degree at most equal to N. We then have, 
for all /: 

(^ajv/)(p) = i^fYip) 
and for all / with degree less than or equal to A — 1: 

(^ajv/)(p)=p-(=^/)(p). 

This relation shows that the operators a^v , ajy after binomial transformation 
turn into the standard creation and annihilation operators [d/dp, p) for a 
restricted set of functions (polynomials of degree at most A). 
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4 Classical dualities of population dynamics 



The scheme developed in the section 2, together with the change of repre- 
sentation discussed in section 3, aUows to recover many of the well-know 
dualities of classical models of population genetics [HHS]. We first consider 
diffusion processes of the Wrigth-Fisher diffusion type and then discrete pro- 
cesses for a finite population of individuals of the Moran type. 

Diffusions and coalescents. 

Consider smooth functions / : [0, 1] — ?> M vanishing at the boundaries and 
1. A diffusion process on [0, 1] has generator of the form 

^ = a{x)^ + (3{x)4- = a{A^)A^ + ^{A^)A , (21) 
dx^ ax 

with A and A^ defined in (jlip . More precisely, we choose 

oo 

a{x) = a^x^ , 

k=l 

oo 

Pix) = 5^/3fcX^ (22) 
where the coefficients a^^Pk satisfy the following 

oo 

"2 = - ^ Ofe, "fe > V/c / 2 , 

fc^2,fe=l 

oo 

/3i = - /3fc, /3fc>0 Vfe/1. (23) 

fc^l,A:=0 

Typical choices are a{x) = x — x"^, f3{x) = (1 — x). By the duality A — a, 
— aT and by theorem 13.11 we find that is dual to 

^f{n) = (^a^a{a^) + al3{a^)^ f{n) 

oo 

= n(n-l) ak{f{n + k-2)-f{n)) 

k^2,k=l 

oo 

+ n I3,{f{n + k-l)-f{n)) (24) 

k^l,k=0 

with duality function D{x,n) = x". By the conditions ()23p on the coeffi- 
cients, this corresponds to a Markov chain on the natural numbers. 
We can then list a few examples. 



14 



1. Wright Fisher neutral diffusion. 

This corresponds to /3 = and —02 = ai = 1 and gives the dual 

^/(n) = {a\a\l-a))))f{n) 

= n(n-l)(/(n-l)-/(n)) , 

which is the well-known Kingman's coalescent block-counting process. 

2. Wright Fisher diffusion with mutation. 

^ = x{l - x)^ + e{l - x)^ = A^l - A^)A^ + ^(1 - A^)A . 
ax^ dx 

This corresponds to ai = —02 = 1, /3o = —f^i = 0. This gives the dual 

^f[n) = [a^{a\l-a^))+6a{l-a))) f{n) 

= n{n - l)(/(n - 1) - f{n)) + en{f{n - 1) - /(n)) , 

which corresponds to Kingman's coalescent with extra rate On to go 
down from n to n — 1, due to mutation. 

3. Wright Fisher diffusion with "negative" selection. 

^ = x{l- x)^ - ax{l - x)4- = A\l - A^){A^ - aA) (25) 
dx^ dx 

with 0" > 0, which corresponds to ai = —02 = 1, /32 = = cr. The 
dual is 

^/(n) = ((a2-cTa)at(l-at))/(n) 

= n(n - l)(/(n - 1) - /(n)) + an{f{n + 1) - /(n)) 

4. Stepping stone model. This is an extension of the Wright-Fisher 
diffusion, modelling subpopulations of which the individuals have two 
types, and which evolve within each subpopulation as in the neutral 
Wright Fisher diffusion, and additionally, after reproduction a frac- 
tion of each subpopulation is exchanged with other subpopulations. 
These subpopulations are indexed by a countable set S. The variables 
Xj G [0, 1], i G S then represent the fraction of type 1 in the i*^ subpop- 
ulation. The generator of this model is defined on smooth functions 
on the set = [0, 1]"^ and given by 
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Here p{i,j), with positive entries outside of the diagonal and with 
^jizsP{h j) = 1) is an irreducible random walk kernel on the set S. 

In terms of the standard creation and annihilation operators = 
Xi,Ai = this generator reads 

^=Y. pi^'jM - - ^i) + E - A^)^^' (27) 

ijGS ieS 

The dual operators act on / : N'^ — )• M via aif{n) = nif{n—ei),a\f{n) = 
f{n + Cj), where 7i € N"^ and Cj G N"^ is the i-th canonical unit 
vector defined via {ei)j = 6ij. By theorem 12.11 item 7, we have 
Ai ai,Al -^^ a] with 

D{n, x) = Y[ ■ 

As a consequence the generator ^ in ()27p is dual to ^ given by 

^= ^pii,j)ia.-aj)ia]-al) + ^aUi^-al) (28) 
i,jes ieS 

or equivalently the generator ^ in ()26p is dual to ^ given by 
■^f{n) = ^p{i,j)ni{f{n-ei + ej) - f{n)) 

+ ^P{j,'i')iT-j{f{n + ei-ej)- f{n)) 
ijes 

+ ^n,(ni-l)(/(7i,-e,)-/(n)) , (29) 

which is the generator of a Markov process on N'^ with transitions n 
n — ei + ej (resp. n ^ n — ej + ej) at rate nip{i,j) (resp. njp{j, i)) and 
n — 7> n— ej at rate ni{ni — 1). The first type of transitions are of random 
walk type and correspond to the exchange of subpopulations, whereas 
the second type are the transitions corresponding to the Kingmans' 
coalescent in each subpopulation. 

Finite-size populations ^ i5j and coalescents. 

As a final example, we illustrate the use of the discrete creation and an- 
nihilation operators a^]^,aN, corresponding to population models with A'" 
individuals in the discrete Moran model. This is the discrete analogue of 
the neutral Wright-Fisher diffusion 

^Nf{k) = — - (^1 - - j (/(fc + 1) + f{k - 1) - 2f{k)) . (30) 
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In terms of the discrete creation and annihilation operators aAr,a]y defined 
in (jlSp . this generator reads 

= ajv(l - a]y)a^ . (31) 

By theorem 12.11 and proposition 13.21 we therefore obtain immediately that 
this generator is dual to the generator of the Kingman's coalescent with 
duality function p^ . 



5 SU{1,1) algebra and corresponding dualities 

In this section we show new dualities for models of population dynamics, 
using dualities between well-chosen differential operators and discrete oper- 
ators. These operators have been used in the context of particle systems 
and models of heat conduction [6]. Interpreted here in terms of population 
models, they yield in that context new dualities. 

We start with the following two families (labeled by m) of infinite dimen- 
sional representations of the algebra SU (1, 1). The first family of operators 
act on functions / : [0, oo) M, whereas the second family acts on functions 
/ : N ^ M. 

= z , 
„ (P m d 

" ^d^^^Tz' 

X» - 4 + (32) 



and 



K^fin) = (y+n)/(n+l) 



K-f{n) = n/(n-l), 

^V(n) = (^+n)/(n). (33) 

The J^f operators satisfy the SU{1, 1) commutation relations 

[^°, = ±J(f^ , [^-, = 2Jf° , (34) 

whereas the K operators satisfy the dual commutation relations (i.e., with 
opposite sign). Therefore, the operators are candidates for a duality relation. 
To find the duality function d : [0, oo) x N — )• M we use 

(d^ 771 d \ 
z-^ + y — j d{z, 0) = K;d{z, 0) = 

which gives as a possible choice (i(z,0) = 1. Then, we can act with J{f~^: 
{Jf.-^rdiz, 0)=z^ = (K+rdiz, 0) = + l) . . . (I + n - l) diz, 77) , 
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and we find 



d{z, n) 



+ 1 



+ n- 1) 



(35) 



One then easily verifies that also d{z, n) = d{z, n) and J^i^d{z, n) = 
K^d{z,n). We can summarize our findings in the following proposition. 

Proposition 5.1. The family of operators given by ()32p and the family 
of operators given by (j33p are dual with duality function given by p5|) . As 
a consequence, every element of the algebra generated by the operators ()32p 
is dual to an element of the algebra generated by ()33p . obtained by replacing 
the operators by their duals and reverting the order of products. 



5.1 Generators from SU{1, 1) raising and lowering operators 

The relevance of the , lies in the fact that some natural generators 
of diffusion processes of population dynamics can be rewritten in terms of 
them. We start now with defining these generators. 

Definition 5.1 ([1], p. 55). The d-types Wright-Fisher model with sym- 
metric parent- independent mutation at rate € M is a diffusion process on 
the simplex Yli=i = 1 defined by the generator 



■t=l * l<i<j<d-l ■' 

'^'^ dg{x) 



1=1 



Definition 5.2. The Brownian Energy process with parameter m € M 
on the complete graph with d vertices (BEP(m)) is a diffusion on Wi^ with 
generator 



^.f-""'/(.) ^ i E yy. (I- - 4) m 

l<i<j<d ^ ^ 



T ^ 

l<i<jd 



Proposition 5.2. The Brownian Energy process with parameter m G M 
on the complete graph with d vertices and with initial condition X^^^^ Xj = 
1 coincides with the d-types Wright- Fisher model with symmetric parent- 
independent mutation at rate 6 = ^{d — 1). 
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Proof. The statement of the proposition is a consequence of the fact the 
BEP evolution conserves the quantity xi + . . . + x^. Consider the initial 
condition Yli=i = 1 ^^id define the function (p : M.'^~^ — > M.^ such that 

d-l 

(Xl,. . . ,Xd„i) = X^ (/)(x) = (Xi, . . .,Xd-l, 1 - ^Xj) . 

i=i 

Then g{x) = f{cj){x)) and, for all i = 1, . . . , d — 1, using the chain rule gives 

dg{x) ^ df{(t>{x)) df{cl){x)) 
dxi dyi dyd 

A computation shows that 

^^^PM f(^xi, . . .,Xd-i,Xd) = ^^rn^^d-i) 9{xi, ■ ■ -^Xd-l) . 

□ 



Definition 5.3. The d-types Moran model with population size N and with 
symmetric parent-independent mutation at rate 9 has generator 

1 
2 



- ^ /ci f fcj + ^— Y j (5f(A; + ei - Cj) +fi((/c - ei + Cj) - 25r(A;)) 

1 "^^^ ( '^^^ 29 \ 

+ 2 ^ ^ " ^ ^ ^^^^ ^ ^ " ""'^ ~ ^^^^^^ ■ 

i=i \ j=i J 

(38) 

Definition 5.4. The Symmetric Inclusion process with parameter m E M 
on the complete graph with d vertices (SlP(m)) is a Markov process on N!|_ 
with generator 

^SlPim) ^^^^ ^ 

I E kJkj + '^){fik + e,-e,) + f{k-e, + e,)-2f{k)). 



2 

l<i<j<d 



(39) 



Proposition 5.3. The generator of the Symmetric Inclusion process with 
parameter m € M on the complete graph with d vertices and with initial con- 
dition X^^^]^ Hi = N coincides with the generator of the d-types Moran model 
with population size N and with symmetric parent-independent mutation at 
rate 9 = f (d- 1). 
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Proof. One verifies that 



SIP{m) 



f{ki,..., kd-l, kd) — -2jvf(i,''^(d-l) 5(^1; • • • 5 kd-l) 



with 



□ 



g{ki,...,kd^i) = f{ki 



d-l 

,kd^i,N -'^kj] 



We can now state our duahty result. 



Theorem 5.1. In the presence of symmetric parent-independent mutation 
at rate 6, the d-types Wright-Fisher diffusion process with generator < IA'6]) 
and the d-types Moran model with N individuals and with generator i38\) 
are dual with duality function 



DNix,k) = 11 



with 



d-l 



d-l 

kd = N-Y,kj . 



(40) 



Proof. The statement of the theorem is a consequence of the duality 
between BEP(m) and SlP(m), which we now recall. We consider the two 
families of operators representing the SU{1, 1) and dual SU{1, 1) commuta- 
tion relations, now rewritten in d coordinates: 



-j^rO _ ™. _8_ I m 
'^m,i — -^i dxi 4 



(41) 



and the corresponding discrete operators 

( K+ ,/(fc,) = (fc, + f - 1) fih - 1) 

K- i/(fe.) = {h + i)f{h + 1) 

I = {h + f ) f{ki) . 

The generator of the BEP(m) then reads 

l<i<j<d ^ 



(42) 



m 



(43) 
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By proposition 15. H combined with theorem I2.H we find that this operator 
is dual to the operator 

L„ = ^ (K.^^u, + K„,K+^^. - 2K^,K5;,,. + ^] , (44) 

l<i<j<d ^ ^ 

This operator is exactly the generator of the SlP(m). The duality function 
is given by, using once more theorem I2.H item 7: 

d 

DN{x,k) = ^d{xi,ki) . 

i=l 

where d{z,k) is given in ()35p . The multiplicative constant r(m/2) in ()35p 
can be dropped, and the result of the theorem thus follows from combining 
the duality between BEP(m) and SlP(m) with proposition 15.21 and proposi- 
tion [531 □ 



5.2 Duality between rf-types Wright-Fisher diffusion and d- 
types Moran model 

We can now also let m — ?• 0, or correspondingly 9 to obtain a duality 
result between the neutral Wright-Fisher diffusion and the standard Moran 
model. 

Let d > 2 be an integer denoting the number of types (or alleles) in a 
population. 

Definition 5.5. The d-types Wright-Fisher model is a diffusion process 
on the simplex XliLi = \ defined by the generator 

i=l * l<i<j<d-l * ^ 

Definition 5.6. The Brownian Energy process with m = on the complete 
graph with d vertices is a diffusion on W^^ given by the generator 

4 _E^^« (I- (40) 

Proposition 5.4. The generator of the Brownian Energy process with 
m = on the complete graph with d vertices and with initial condition 
Yli=i Xi = I does coincide with the generator of the d-types Wright- Fisher 
diffusion, i.e. 

-^/^^^°V(a;i, • • .,Xd-i,Xd) = ^d^^g{xi, . . .,Xd-i) 
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with 

d-l 

g{xi, . . .,Xd-i) = f{xi,.. ■,Xd-i, 1 - ^^Xj) 

3=1 

Proof. Similar to the proof of proposition 15.21 □ 



Definition 5.7. In the d-types Moran model with population size N a 
pair of individuals of types i and j are sampled uniformly at random, one 
dies with probability 1/2 and the other reproduces. Therefore, denoting by 
k = {ki, . . . ,kd-i), where ki is the number of individuals of type i , the process 
has generator 



■^N,Tg{k) = \ k,kjigik + e,-ej)+g{k-ei + e,)-2g{k)) 

l<i<j<d~l 

+ i^^h iV-^fc,- {g{k + ci) + g{k - Ci) - 2g{k)) . 



i=i \ j=i 



(47) 



Definition 5.8. The Symmetric Inclusion process with m = on the 
complete graph with d vertices is a Markov process on N!^ with generator 

= I Yl kikj{f{k + e,-ej) + fik-e, + ej)-2fik)). 



2 

l<i<j<d 



(48) 



Proposition 5.5. The generator of the Symmetric Inclusion process with 
m = on the complete graph with d vertices and with initial condition 
Yli=i "'j = does coincide with the generator of the d-types Moran model 
with population size N, i.e. 

^(°V(fci, ■ ■ ■ , kd-iM = ^If/giki,. . . , kd-i) 

with 

d-l 

g{ki,...,kd^i) = f{ki,...,kd-.i,N -^kj) . 

i=i 

Proof. Similarly to the proof of proposition 15.31 the result follows from 
the conservation law, namely the fact that the SIP evolution conserves the 
total number of particles ki + . . . + kd- □ 

In the duality result of theorem 15.11 we cannot directly substitute m = 
because there would be problems when some ki = 0. To state a duality 
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result for 6 = 0, i.e., between the Wright Fisher diffusion and the Moran 
model without mutation, we start again from the duality between Brownian 
Energy process and Symmetric Inclusion process: 



IEf^"^™M n 7.^^?^^ (49) 



jg57P(m) 




For ^ € N'', denote = {i : d > 1}. Then we can rewrite (gO]) and 

obtain 

(d ^i{t) 
JL (f+i)...(f+?.w-i) 

Notice that in the SIP{0) process, only coalescence jumps take place, hence 
^{iit)) < for all t > 0, ^ G N'^. Now we are in the position to take 

the limit m — ?> and we find 

eBep,„,| n^^^j (51) 

Ef I / («(£) = «({(())) n 



Notice that the Ihs becomes zero as soon as one of the Xi are zero, which 
corresponds to the fact that for all i, Xj = is an absorbing set in the 
diffusion. Corresponding to this, the rhs becomes zero as soon as one of the 
species disappear, i.e., as soon as ^(^) decreases by one unit. We call ([5T]) 
^''duality until absorption with duality function" 

mx)={ n 7#I)t| • (52) 

We can summarize our findings in the following theorem. 

Theorem 5.2. The d-types Wright-Fisher diffusion process with generator 
(T^y and the d-types Moran model with N individuals and with generator 
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I^Tp are dual until absorption on duality function 




(53) 



with 



d-1 



d-1 



'3 ' 



kd = N 



Remark 5.1. We remark that the duality results in subsections \ 5.1\ and 
\5.2\ are of a different nature than the usual dualities between forward pro- 
cess and coalescent. Indeed, we have here duality between two "forward pro- 
cesses" (the Wright-Fisher diffusion and the Moran model), which cannot be 
obtained from "looking backwards in time", the method by which moment- 
dualities with the coalescent are usually obtained. In our framework, the 
dualities with the coalescent correspond to a change of representation in the 
Heisenberg algebra, whereas the dualities between e.g. Wright-Fisher and 
Moran model arise from a change of representation in the SU{1, 1) algebra. 

5.3 Self-duality of the d-types Moran model 

We can push further the SU{1, 1) structure behind the Moran model and 
deduce self-duahty of the process. 

Theorem 5.3. The d-types Moran model with N individuals and with gen- 
erator 123) is self-dual with duality function 



Proof. The result follows from the self-duality property of the SlP(m) 
process [8] and from proposition 15.31 □ 

The limit m — )• 0, or equivalently ^ ^ 0, leads to the self-duality until 
absorption relation 




(54) 



where kd = N - Ya=i h and = N - Yli=i d- 




(55) 



This computation then leads to the following theorem. 
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Theorem 5.4. The d-types Moran model with N individuals and with gen- 
erator j-^^Tp is self-dual until absorption, with duality function 



where = N - YliZl h and = N - Y^i~l Ci- 
5.4 Examples 

We give some examples of concrete computations using the dualities of the 
present section. 

1. Heterozygosity of two- types Wright-Fisher diffusion. This is exactly 
the BEP{0), xt, yt on two sites, with initial condition x,y : x + y = 1. 

^^f^'\x{t)y{t)) = Eff(°)(xy/(ni(t) = l,n2(t) = l)) 

= xyPi,i(ni(t) = 1, n2{t) = 1) = xye'* , 

where denotes the law of the 5/P(0) process initialized with one 
particle per site. 

2. Higher moments of two-types Wright-Fisher diffusion. We use the 
same notation of the previous item and consider for instance x'^y. Fur- 
ther, we notice that if we start the SIP{0) from initial configuration 
(ni, 77-2) = (2, 1), then the only transitions before absorption are of the 
type (2, 1) (1,2) and vice versa, and both transitions occur at rate 
2, whereas from any of these states, the rate to go to the absorbing 
states is also equal to two. Therefore, 

E%^^('Hx\t)y{t)) 

= x2yPf,^(°)((ni(i),n2(t)) = (2, 1)) + xy^ff (°)((ni(t), n2(t)) = (1, 
= ^(:rVl + e-^*)+xy2(l_e-2*)). 



3. Analogue of heterozygosity for d-types Wright Fisher diffusion. Notice 
that for multitype Wright Fisher, there is no simple analogue of the 
Kingman's coalescent, as for the two- types case. This means that we 
have the BEP{0) started from xi, . . . ,Xd 

= XI . . . XdE^J^l^^l^ (/(n,(t) / Vi G {1, . . . , d})) 
= XI . . . x.e-^'^-^)* . 
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4. Analogue of x^y for the multi-type case. 



<S°](^?W^2...X,(t)) 

= XI . . . XdE^^l^^^l^ {I{m{t) / Vi G {1, . . . , d})) 



d 



= E n^i ^?i'f2X = 1' • • • '^^(*) = 2' ■ ■ -^Mt) = 1) 

To compute the latter probability, we remark that starting from the 
configuration (2, 1, ... , 1), the SIP{0) will be absorbed as soon as one 
of the particles on the sites with a single occupation makes a jump, 
which happens at rate (d — l)((i — 2) + (d — 1)2 = d(d — 1). Further, as 
long as absorption did not occur, the site with two particles moves as a 
continuous-time random walk Xf on the complete graph of d vertices, 
moving at rate 2 and starting at site 1. Therefore 

^tli!^Z)iMt) = !,•••, riiit) = 2, . . . , nd{t) = 1) 

= e +-(l-e (1 - - e j. 

As these examples illustrate, computations of moments in the multi-type 
Wright Fisher model reduce to finite dimensional Markov chain computa- 
tions, associated to inclusion walkers on the complete graph until absorption, 
which occurs as soon as a site becomes empty. The same can be done for 
the multi-type Moran model, using its self-duality. 
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